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Abstract 

Background: Sturgeon cultivation is important for botli industry and aquaculture in Cliina. To date, more tlian 17 
species or strains liave been farmed for fillets and caviar production. Crossbreeding among different sturgeon 
species is frequent and the F2 hybrids are fertile. However, large-scale farming can have negative impacts on wild 
populations i.e. escape of exotic sturgeons and must be taken into consideration. Escape of exotic sturgeons can 
cause severe ecological problems, including threatening native sturgeon species once the exotic varieties become 
established or hybridize with native individuals. However, little is known about their genetic resources and variation. 

Methods: Genetic diversity and introgression of seven sturgeon species were analyzed using mitochondrial DNA 
cytochrome oxidase subunit I (CO/) and nine microsatellite markers. This study included 189 individuals from seven 
sturgeon species and 277 individuals from ten lineages of F2 hybrid strains. 

Results: MtDNA CO/ sequences (632 bp long) were generated from 91 individuals across the 17 sturgeon strains 
and produced 23 different haplotypes. Haplotype diversity was high (h = 0.915 ± 0.015) and nucleotide diversity 
was low (tt = 0.03680 ± 0.00153) in the seven sturgeon species and ten interspecific hybrids. Phylogenetic analyses 
resulted in almost identical tree topologies, and different haplotype structures were mainly related with sturgeons 
of different female parents. Analysis of molecular variance revealed that 81.73% of the genetic variance was due to 
matrilineal differences, while 9.40% resulted from strain variation. Pairwise Fst values obtained with POLYSAT 
software, were high among strains and ranged from 0.031 to 0.164. Admixture analysis assigned seven distinct 
groups and ten genotypes of admixed clusters composed of hybrid strains using STRUCTURE when assuming K = 7. 

Conclusions: The interspecific mtDNA gene tree corresponded to the expected taxonomic divisions. These 
relationships were also supported by the results from the microsatellite analysis and contributed to unambiguously 
identify seven sturgeon species and ten F2 hybrid strains from sturgeon farms in China. Moreover, we found that 
introgressive hybridization is pervasive, exists in both purebred and hybrid sturgeons, and may reflect widespread 
mismanagement in sturgeon breeding in China. 



Background research has been performed to evaluate crossbreeding 

Hybridization is probably an inevitable process during in farm animals [10-12]. Crossbreeding is fast and effect- 

speciation [1], in which offspring inherit restructured ive, creates heterosis in the future generations [13,14] 

parental genes obtained by mating individuals of differ- and can contribute to genetic improvement [15]. 

ent genotypes. Hybridization is common both in plants Sturgeon is one of the most ancient fish in the world 

and animals [2-4], and is largely exploited to improve and belongs to the order Acipenseriformes that contains 

various species [5-7], especially plants [8,9]. Through 27 species divided into two families i.e. Acipenseridae 

hybridization, desirable traits can be combined, leading (sturgeon, 25 species) and Polydontidae (paddlefish, two 

to more competitive descendants and heterosis. Much species). Most sturgeon species are near extinction [16], 

and have been listed in the Appendices to CITES (Con- 

— \ \ . — , vention on International Trade in Endangered Species of 
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geons are considered worldwide as excellent candidates 
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for aquaculture due to their high commercial value (cav- 
iar production and meat). Conflicting issues between the 
drastic reduction in natural populations and the huge 
profits in the business of sturgeon farming are the main 
driving force promoting the development of sturgeon 
aquaculture in China [17]. Recently, with the develop- 
ment of sturgeon aquaculture and the increase in 
cultured broodstock, China has become the largest stur- 
geon aquaculture country in the world, and hybrid 
sturgeons are widely bred [17]. The most recent survey 
has recorded more than 17 sturgeon strains (including 
both purebreds and hybrids) aquacultured in China, with 
two species i.e. Acipenser baerii and A. schrenckii, and 
three hybrid sturgeon strains i.e. A, baerii x A, 
schrenckiiy A, schrenckii x Huso dauricus, and A, baerii x A, 
gueldenstaedti dominating sturgeon farming [17]. Among 
the hybrid sturgeon lines, some grow slowly and are of 
low quality value, but it is difficult to identify hybrid 
strains of sturgeon, especially at the early stages of 
growth i.e. fry and sub-adults. In addition, little is 
known on the genetics and polymorphism of stur- 
geons. In recent years, progress in artificial propaga- 
tion technology has expedited the development of 
sturgeon aquaculture. However, when selective breed- 
ing is performed, the genetic relationships between 
candidate parents are unknown, resulting in random 
hybridization of different sturgeons. To prevent both 
inbreeding and degradation of germplasm resources, it 
is necessary to analyze the genetic background of the 
candidate breeders. Mitochondrial DNA (mtDNA) 
and microsatellite genetic markers have been used 



successfully in the studies of fish for introgressive 
hybridization [18-21]. In this study, we used nine 
microsatellites and the mitochondrial DNA (mtDNA) 
gene cytochrome oxidase subunit I (COI) to examine 
the genetic diversity, introgression, and differences in 
genetic background of 17 farmed sturgeon strains, in- 
cluding seven species and ten hybrids. 

Methods 

Animals and DNA extraction 

From 2008 to 2011, seven sturgeon species and ten 
hybrids were obtained from the Engineering and Tech- 
nology Center of Sturgeon Breeding and Cultivation of 
Chinese Academy of Fishery Science (Beijing, China), 
the Heilongjiang Fisheries Research Institute, the Beijing 
Fisheries Research Institute, and the Hangzhou Qiandaohu 
Xunlong Sci-Tech Co., Ltd. (Table 1). The historical 
background of all the sturgeons collected was known, and 
each species was identified according to anatomical charac- 
teristics. Fish fins were stored in 95% ethanol. Genomic 
DNA was prepared using a DNeasy tissue kit (Qiagen). 

Mitochondrial DNA gene COI amplification and 
microsatellite genotyping 

The mtDNA gene COI (632 bp) was amplified for five to 
eight individuals of each sturgeon strain (total n = 91) using 
published primers that were designed specifically for fish 
i.e. (FishFl [5'-TCAACCAACCACAAAGACATTGGCA 
C-3'] and FishRl [5'-TAGACTTCTGGG TGGCCA 
AAGAATCA-3']) [22]. PGR (polymerase chain reaction) 
amplifications were carried out in 25 [iL reaction volumes 



Table 1 Sturgeon species and code, number of Individuals and line type used in the experiment 



Species 


Species code 


Number of individuals 


Type 


A. schrenckii 


S 


30 


Purebred 


A. boerii 


X 


30 


Purebred 


A. gueldenstaedti 


E 


34 


Purebred 


H. dauricus 


H 


26 


Purebred 


A. rutlienus 


Xi 


23 


Purebred 


A. sinensis 


Z 


25 


Purebred 


A. stellatus 


G 


21 


Purebred 


A. boerii'^ x H. douricusS 


XH 


31 


Hybrid 


A. boerii^ x A. scinrencl<iiS 


XS 


32 


Hybrid 


A. schrencliii^ x A. boeriiS 


SX 


30 


Hybrid 


A. boerii^ x A. gueldenstaedti S 


XE 


23 


Hybrid 


A. gueldenstaedti^ x A. baeriiS 


EX 


40 


Hybrid 


A. gueldenstaedti^ x H. dauricusS 


EH 


30 


Hybrid 


H. dauricus'^ x A. schrenckii^ 


HS 


30 


Hybrid 


A. schrenckii'^ x H. dauricusS 


SH 


26 


Hybrid 


A. ruthenus'^ x H. dauricusS 


XiH 


30 


Hybrid 


A. sinensis^ x A. schrenckii S 


ZS 


5 


Hybrid 
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containing approximately 20 ng of template DNA, 1 (iL of 
each primer, 9.5 (iL of ddH20, and 12.5 (iL of 2x Premix 
Taq DNA polymerase (Takara). Amplification conditions 
were as follows: 95°C for 2 min; 35 cycles at 94°C for 30 s, 
57°C for 30 s, and 72°C for 1 min; and a final extension at 
72°C for 10 min. Positive (with DNA template) and nega- 
tive (with water) controls were used to check PGR perform- 
ance and contamination. The PGR products were purified 
using the PGR purification kit (Shanghai Bio-Tec, Ltd) and 
sequenced with the ABI PRISM BigDye Terminator Ready 
Reaction kit (Applied Biosystems) and run on an ABI 377 
genetic analyzer according to the manufacturer s protocol. 
To avoid errors in sequencing, PGR amplifications of all 
samples were sequenced on both strands. 

Fifteen microsatellite loci were selected before the start 
of the experiment, the amplified products were sepa- 
rated by electrophoresis on 10% (29/1 Acrylamide- 
bisacrylamide) 0.50 mm thick denaturing polyacrylamide 
gels (120 V, 10 h), and DNA bands were visualized using 
silver staining. The size of individual alleles was deter- 
mined using a 50 bp DNA size standard (Tiangen) and 
allele sizes were estimated by Gelpro 3.2. Nine of 15 
most effective loci were selected for further analyses 
after the preliminary analysis [See Additional file 1: 
Table SI] [23,24]. The corresponding nine primer pairs 
were fluorescently labeled with one of the following dyes 
FAM, HEX, or TAMRA. PGR was performed under the 
following conditions: preheating at 94°G for 4 min; 
denaturing at 94°G for 40 s, annealing at 50-58°G for 30 
s, and elongation at 72°G for 40 s for 35 cycles; exten- 
sion at 72°G for 10 min. Amplifications were carried out 
in 15 (iL reaction volumes containing 10 ng of template 
DNA, 0.5 (iL of each primer, 6 (iL of ddH20, and 7.5 (iL 
of 2x Premix Taq DNA polymerase (Takara). All PGR 
were carried out in a PTG-100 thermal cycler (MJ Re- 
search). Allele sizing was carried out by automated 
fluorescent scanning detection in an ABI 377XL DNA 
sequencer (Applied Biosystems) using ROX500 as in- 
ternal lane size standard, and the software GENE SG AN 
and GENOTYPER (Applied Biosystems). 

Sequence analysis 

Nucleotide sequences were aligned using the GLUSTALX 
software (1.83) and checked by eye. Variable nucleotide 
sites (including transversions (tv), transitions (ts) and inser- 
tions/deletions) were calculated using MEGA 2.1 software. 
Haplotypes were estimated using DNASP, and the haplo- 
type network was drawn using NETWORK software. For 
the phylogenetic analysis, Modeltest 3.06 was run to deter- 
mine the appropriate model of molecular evolution in a 
likelihood ratio test framework. We then performed 
maximum-parsimony (MP) analysis using the programs 
PAUP 4.0 and MRBAYES 3.0, with Polypterus bichir bichir 
as out-group. Gaps were treated as missing in the 



parsimony analyses. Bootstrap analyses were performed 
with 6000 replicates and 1000 fiall heuristic replicates for 
maximum parsimony. For Bayesian phylogenetic inference, 
four Markov chain Monte Garlo (MGMG) simulations were 
run for 1 500 000 generations, sampling every 1000 genera- 
tions. The initial 5% of trees were discarded as burn-in, and 
finally, a 50% majority rule consensus tree was constructed. 

Specific nucleotide sites were analyzed in the COI 
sequences, and haplotype distributions, nucleotide diver- 
sities (71) within and among strains were estimated using 
DNASP. Hierarchical analysis of molecular variance 
(AMOVA) was performed to compare the levels of genetic 
diversity within and among several possible strain groups 
using Arlequin 3.11 with 1023 permutations. 

Microsatellite analysis 

Individual polyploidy genotypes were scored from 
microsatellite banding patterns in the electropherograms 
according to the Microsatellite DNA Allele Gounting- 
Peak Ratios (MAG-PR) method of Esselink et al.(2004) 
[25]. The STRUGTURE version 2.3.2 program [26] was 
used to infer strain composition based on microsatellite 
data. We used the Bayesian clustering approach in 
STRUGTURE 2.0 to identify the most likely number of 
clusters (K) as well as to assign individuals to these clus- 
ters. We performed five replicate runs (burn-in period of 
1 000 000 steps and 120 000 MGMG iterations) at each 
value of K from 2 to 7. Individuals were assigned on the 
basis of their membership coefficient and factor correl- 
ation analysis (EGA) based on the multilocus genotypes 
was carried out using the GENETIX program to separate 
the strains and identify any intermediate genotypes 
resulting from admixture of strains [27]. 

We used the POLYSAT software to calculate the num- 
ber of alleles, allele frequency, and Shannon- Wiener 
index. The Simpson index was used to characterize the 
levels of genetic diversity in each sturgeon strain [28]. 

Results 

Mitochondrial DNA 

MtDNA COI sequences (632 bp) were generated from 17 
strains represented by 91 individuals. One hundred and 
six variable nucleotide sites were found, including 14 
transversions (tv), 82 transitions (ts) and 10 insertions/dele- 
tions, and 23 haplotypes were defined [GenBank accession 
numbers: KG578823-KG578845]. Our results revealed high 
haplotype diversity (h = 0.915 ± 0.015) and low nucleotide 
diversity (tt = 0.03680 ± 0.00153) in the seven sturgeon spe- 
cies and ten interspecific hybrids. The seven pure sturgeon 
species harbored unique nucleotide sites: A. ruthenus (Xi) 
had the greatest number of specific nucleotide sites (i.e. 14), 
while A. baerii (X) and A. gueldenstaedti (E) had the lowest 
number (i.e. 1) [See Additional file 2: Table S2]. H09 was the 
most common haplotype shared by 19 individuals and A. 
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ruthenus^ x H. dauricusS (XiH) had the biggest number of 
haplotypes (Table 2). 

Using MODELTEST, the HKY+G model was built 
using the best fitting distance estimator. It showed a 
gamma distribution shape parameter of 0.2655, a transi- 
tion/transversion (Ti/Tv) ratio of 4.0901, and base fre- 
quencies of A = 0.2650, C = 0.2825, G = 0.1780 and T = 
0.2744. Parsimony analyses were performed under equal 
weight (Ti/Vi = 1) and unequal weight (Ti/Vi = 4.0901) 
sets. All phylogenetic analyses resulted in almost identi- 
cal tree topologies. With the exception of all individuals 
from A. baerii^ x A, gueldenstaedtiS (XE) and one indi- 
vidual from A, ruthenus^ x H, dauricusS (XiH), most 
of the tree topologies were grouped together with the 
same matrilineal sturgeon [See Additional file 3: Figure 
SI]. This observation was supported by haplotype net- 
work analysis (Figure 1). 

AMOVA analysis revealed that 81.73% of the genetic 
variance was due to matrilineal differences and 9.40% to 
strain variation, while only 8.86% of the genetic variance 



originated from strains within groups with the same 
matrilineal sturgeon (Table 3). 

Microsatellites 

A total of 466 individual sturgeons (Table 1) were geno- 
typed using nine nuclear microsatellite markers [See 
Additional file 4: Table S3]. All tested microsatellite loci 
were polymorphic in all farmed sturgeons. Microsatellite 
LS68 had the biggest number of alleles per individual, while 
LS19 had the smallest. A, sinensis^ x A, schrenckiiS (ZS) 
tended to have more alleles per individual per locus, while 
A, baerii^ x H, dauricusS 0^) had the fewest [See 
Additional file 5: Table S4]. 

Using POLYSAT software, we showed that the total 
number of alleles observed at each locus ranged from 7 
to 15, with a total of 386 alleles detected across the nine 
loci. The total number of alleles found within strains 
ranged from 42 {A, sinensis^ x A, schrenckiiS) to 117 
(//. dauricus^ x A, schrenckiiS) y and the number of 



Table 2 Summary of mtDNA COI region haplotype distribution among the 17 sturgeon species 

Species XH X XS SX S SH XE EX E EH H H XiH Xi ZS Z G Total 

Haplotype 



HOI 1 1 

H02 2 2 4 

H03 2 3 5 10 

H04 5 5 

H05 1 1 

H06 1 1 

H07 1 1 

H08 2 2 

H09 6 4 5 4 19 

H10 1 1 

H11 1 1 

H12 1 1 

H13 4 4 8 

H14 1 1 

H15 5 16 

H16 5 5 

H17 1 1 

H18 1 1 

H19 5 5 

H20 5 5 10 

H21 1 1 

H22 4 4 

H23 2 2 

Total 5 55555 6555 558 5566 91 



Code for species is presented in Table 1. 
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0 A baerii 

^ A gueldenstaedti 
0 A ruthenus 
0 A schrenckii 
Q A sinensis 
Q H. dauricus 
0 A stellatus 



Figure 1 Statistical parsimony networic based on tiie mtDNA gene CO/ iiaplotypes. Each circle represents a single haplotype; circle size is 
scaled by haplotype frequency; the same colors indicate the same matrilineal sturgeon; 1: /\. boerii, A. boerii^ x H. dauricusS and A. baerii'^ x A. 
schrenckiiS; 2: A. gueldenstaedti, A. gueldenstaedti'^ x A. baeriiS, A. gueldenstaedti'^ x H. dauricusS and A. baerii'^ x A. gueldenstaedtiS; 3: A. 
ruthenus and A. ruthenus^ x H. dauricusS', 4: A. schrenckii, A. schrenckii^ x A. baeriiS, A. schrenckii^ x H. dauricusS and A. ruthenus^ x H. 
dauricusS', 5: A. sinensis and /\. sinensis^ x /\. schrenckiiS', 6: /-/. dauricus and /-/. dauricus^ x /\. schrenckiiS; 7: /\. stellatus. 



alleles per microsatellite ranged from 25 (LS54) to 71 
(SPL106) [See Additional file 4: Table S3]. Shannon- 
Wiener and Simpson indices were used as measures of 
allelic diversity and dominance in the population, re- 
spectively. Based on the microsatellite analysis, the 
Shannon- Wiener and Simpson indices of sturgeon alleles 
ranged from 1.561 to 3.404 and from 0 to 0.092, respect- 
ively (Figure 2). 

Factor correlation analysis (FCA) based on microsatel- 
lite genotypes revealed a clear segregation among 



purebred sturgeon species and hybrids. FCA highlighted 
the differences among the seven analyzed purebred spe- 
cies. The most informative is axis 1 (23.18% of the total 
genetic variation), which separates the seven species and 
distinguishes A, gueldenstaedti (E) and A. baerii (X) 
from the five other species. Axis 2, which is slightly less 
informative (15.54%), separates A, ruthenus (Xi), A, 
sinensis (Z), A, schrenckii (S) and A, stellatus (G). Finally, 
axis 3 mainly separates A. gueldenstaedti (E) and H, 
dauricus (H) with 12.07% of the total genetic variation 



Table 3 AMOVA results for sturgeon CO/ gene estimated using F-statistics 



Source of variation d.f. Sum of squares 

Among groups 6 910.209 

Among populations within groups 10 76.851 

Within populations 74 95.875 

Total 90 1082.934 
* P < 0.01 ; d.f. = degree of freedom. 



Variance components 

1 1 .26227 
1.22116 
1.29561 
13.77904 



% of variation 

81.73 

8.86 

9.40 



F-statistics (P-value) 

FST:0.90597* [P = 0.00000) 
FSC:0.48521*(P = 0.00293) 
FCT:0.81735* = 0.00000) 
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Figure 2 The Shannon-Wiener and the Simpson indices of 17 sturgeon strains. 



(Figure 3). The ten hybrids, even with their typical 
mtDNA COI haplotypes of female descent, are identified 
as hybrids, but occupy an intermediate position among 
the parents according to the multidimensional analysis of 
microsatellite genotyping data [See Additional file 6: Figure 
S2, Additional file 7: Figure S4, Additional file 8: Figure S5, 
Additional file 9: Figure S6, Additional file 10: Figure S7, 
Additional file 11: Figure S8 and Additional file 12: 
Figure S9). 

To investigate the existence of different genetic 
clusters in sturgeon breeds, the admixture model 
implemented by the STRUCTURE software was used 
to analyze separately the microsatellite genotyping 
data. The assignment test, which was conducted with 
STRUCTURE with K ranging from 2 to 7 and was 
based on nine microsatellites, separated the samples 
into seven distinct groups, congruent with the seven 
sturgeon species, and ten mixed-cluster groups 
representing the ten interspecific hybrid strains [See 
Additional file 13: Figure S3]. Various degrees of intro- 
gressive hybridization were also observed in all seven 
purebred sturgeons with H, dauricus showing the 
highest level. 



Discussion 

Molecular identification 

DNA barcoding is widely used for species determination 
because sequence divergences are generally much lower 
among individuals of a given species than among closely 
related species [29-32]. However, hybridization among 
species can cause taxonomic uncertainty. Since mtDNA 
is maternally inherited, any hybrid or subsequent gener- 
ation will only have maternal mtDNA [31]. 

Jenneckens et al. (2000) reported that genetic contam- 
ination of A, gueldenstaedti with A, baerii or A, baerii 
hybrids occurred in the Volga River in Russia. Crosses 
and backcrosses of these species with native A, 
gueldenstaedti led to the loss of the morphological diag- 
nostic A, baerii features [33]. In our study, we also used 
nine microsatellite loci from nuclear DNA and mtDNA 
COI for sturgeon identification. Our data showed that 
mtDNA can identify purebred sturgeon species and the 
maternal origin of hybrids (Figure 1) and [See Additional 
file 3: Figure SI] since each line harbors unique sites 
[See Additional file 2: Table S2]. The interspecific 
mtDNA gene tree constructed agreed with the expected 
taxonomic divisions, except for all individuals of A, 
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baerii^ x A, gueldenstaedtiS (XE) and one individual of 
A, ruthenus^ x //. dauricusS (XiH). A, baerii^ x A, 
gueldenstaedtiS (XE) was grouped together with strains 
bred from the female parent A, gueldenstaedti, perhaps 
because of sex-specific or nonreciprocal crossing. 
Moreover, A, baerii may be a hybrid and may have 
been backcrossed many times with A, gueldenstaedti. 
Hybridization may be involved in particular gender 
combinations in each species due to the maternal in- 
heritance of mtDNA [34]. One individual of A, 
ruthenus^ x H, dauricusS (XiH), identified within the 
wrong cluster, may be due to erroneous identification 
of the original specimen. 

In general, matrocliny and patrocliny with different 
combinations of dominant and recessive alleles of the 
genes determine the diagnostic characters. These clearly 
demonstrate how the hybrid individuals acquire varying 
degrees of similarity and differences from the parental 
species [35]. Taxonomic confusion as a result of inter- 
specific hybridization does not seem to be a major issue 
and can be solved with genetic markers of matroclinous 
and patroclinous inheritance [36]. Mitochondrial DNA 



sequences are frequently transferred to the nucleus, 
giving rise to the so-called nuclear mtDNA sequences 
or NUMT (nuclear mitochondrial DNA) because of 
hybridization [37]. NUMT is not helpful for species 
identification because of its potential pseudogene status 
[31]. NUMT are common in plants and animals, but few 
NUMT are found in fish [38,39]. In our results, since no 
evidence of significant linkage disequilibrium (LD) was 
found among any loci, we can avoid this negative factor 
of NUMT for random distribution of these loci. 
Microsatellite results can distinguish reciprocal crosses 
among different sturgeons [See Additional file 6: Figure 
S2, Additional file 7: Figure S4, Additional file 8: Figure 
S5, Additional file 9: Figure S6, Additional file 10: 
Figure S7, Additional file 11: Figure S8 and Additional 
file 12: Figure S9), and a few samples were assigned to 
the wrong' congeneric species, which may represent 
introgressive hybridization. 

Introgression 

Microsatellite analysis reveals population structure and 
helps to explain past introgression of DNA among 
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sturgeon species. Strong signatures are widely docu- 
mented when comparing introgression of microsatellites 
with that of mtDNA [36,40]. We found that introgres- 
sive hybridization is pervasive in either purebred or hy- 
brid farmed sturgeons, which may reflect widespread 
mismanagement of sturgeon breeding in China. STRUC- 
TURE results showed that introgressive hybridization in 
domestic varieties (A. schrenckii, H, dauricus, and A. 
sinensis) is less serious than in varieties introduced to 
China from abroad {A, baerii, A, gueldenstaedti, A. 
ruthenusy and A, stellatus), H, dauricus is the most obvi- 
ous example of strain impurity with a high level of gen- 
etic information originating from another sturgeon [See 
Additional file 13: Figure S3], which may be due to the 
fact that this species was frequently used as a parent in 
crosses. Our results also showed that A, stellatus was 
almost assimilated by A, baerii and that A. ruthenus had 
introgressed into the hybrid stains of A, baerii as a 
parent, especially in A. baerii^ x A. schrenckiiS and 
A, gueldenstaedti^ x A. baeriiS- Such sturgeon 
hybridization has also been reported between the en- 
dangered A. ruthenus species and the exotic A. baerii 
species in the Danube River and between A. baerii and 
A. ruthenus in the Irtysh River of Sinkiang [41,42]. 
Campton reported that introgression can cause a more 
abundant species to genetically assimUate a rare spe- 
cies [43]. 

In China, it is easy for farmed sturgeons to escape 
to natural water systems where wild sturgeons live, and 
there are no behavioral or physiological barriers 
preventing interbreeding between the cultured and wild 
populations [44]. A, sinensis is distributed in the Yangtze 
River Basin, but recently, different types of sturgeons 
have been cultured along the tributaries, and escaped in- 
dividuals have crossed with A, sinensis ^ leading to its 
farraginous genotype. 

Genetic diversity 

MtDNA results revealed high gene diversity and low 
nucleotide diversity in these farmed sturgeons. Results 
from AM OVA analysis for the 17 strains of sturgeon 
showed that most of the genetic variance was distributed 
among the different matrilineal strains. These conclusions 
were supported by our microsatellite results. Since 
Shannon- Wiener and Simpson indices measure allelic gene 
diversity and dominance in the population, respectively, a 
high and homogeneous Shannon- Wiener index may indi- 
cate that hybrid regions could comprise a large number of 
multigenerational backcrossed hybrids that are indistin- 
guishable from the parental species [45], and that multi- 
generational backcrossed hybrids resemble the parental 
species with which they share most of their alleles [46]. A, 
ruthenus^ x H, dauricusS (XiH) showed the highest 
Simpson index, which may be due to the fact that 



introgression of a few loci may promote adaptive 
divergence [1]. Similarly, high levels of genetic diversity 
have been reported in lake sturgeons across the species' 
range [46]. The low levels of genetic diversity in A, 
sinensis^ x A, schrenckiiS (ZS) suggest that A. sinensis^ x 
A. schrenckiiS (ZS) was founded by a small number of 
individuals. 

Genetic diversity of farm animal is often exploited to 
meet current production needs, to allow sustained genetic 
improvement, and to facilitate rapid adaptation to chan- 
ging breeding objectives [47]. Allele introgression is crucial 
to livestock genetics, while crossbreeding has always been 
a staple of breeding programs. Hybridization and polyploid 
formation will continue to generate species diversity . Our 
results show that farmed sturgeons hybridize with each 
other to varying degrees and with high genetic diversity. 
The conservation aquaculture programs are designed to 
minimize the genetic impacts on wild populations caused 
by introgressive farmed populations. 

Conclusions 

We combined a morphological study with DNA barcoding 
and microsatellite markers to show that we can unambigu- 
ously identify seven sturgeon species and ten F2 hybrid 
strains from sturgeon farms in China. We found that intro- 
gressive hybridization is pervasive and exists in both pure- 
bred and hybrid farmed sturgeons, which may reflect 
widespread mismanagement of sturgeon breeding in China. 
We also found that farmed sturgeons hybridize with each 
other to varying degrees and with high genetic diversity. 

Additional files 



Additional file 1: Table SI. Characteristics of sturgeon microsatellite 
markers in this study. The data provided represent characteristics of the 
nine most effective microsatellite loci which were selected for analyses in 
this work. 

Additional file 2: Table 52. Specific nucleotide sites observed in the 
mtDNA COI sequences of seven purebred sturgeon species. The data 
provided represent specific nucleotide sites analyzed by MEGA, all the 
sites were observed in mtDNA COI sequences of seven purebred 
sturgeon species. 

Additional file 3: Figure SI. Phylogenetic tree obtained by Bayesian 
inference from the analysis of COI sequences with Polypterus bichir bichir 
as out-group. The figure represents the phylogenetic tree infered from 
the analysis of COI sequences with Polypterus bichir bicliir as out-group. 
Different colors represent different matrilineal. 

Additional file 4: Table S3. Total number of alleles detected for the 
nine microsatellite loci selected for the analysis of 17 sturgeon strains. 
The data provides the number of alleles detected for the nine 
microsatellite loci by POLYSAT. 

Additional file 5: Table S4. Average number of alleles per individual 
per locus detected in 17 sturgeons. The data provides the average 
number of alleles per individual per locus detected in 17 sturgeon strains. 
N is the total number of individuals for each sturgeon strain. 

Additional file 6: Figure S2. Factorial Correspondence Analysis (FCA) 
based on nine microsatellite loci in A. baerii, A. gueldenstaedti and their 



Zhang et at. Genetics Selection Evolution 2013, 45:21 
http://www.gsejournal.0rg/content/45/l/2i 



Page 9 of 10 



hybrids. The figure shows the distir^ction betweer^ A. baerii, A. gueldenstaedti 
and A. baerii^ x A. gueldenstaedti^ ar^d A. gueldenstaedti^ x A. baerii^. 

Additional file 7: Figure S4. Factorial Corresponder^ce Analysis (FCA) 
based on nine microsatellite loci in A. schrenckii, H. dauricus and their 
hybrids. The figure shows the distinction between A. schrenckii, H. 
dauricus and A. schrenckii'^ x H. dauricusS and H. dauricus^ x A. 
schrenckii S- 

Additional file 8: Figure S5. Factorial Correspondence Analysis (FCA) 
based on nine microsatellite loci in A. baerii, A. schrenckii and their 
hybrids. The figure shows the distinction between A. baerii, A. schrenckii 
and A. baerii^ x A. schrenckii^ and A. schrenckii^ x A. baerii^. 

Additional file 9: Figure S6. Factorial Correspondence Analysis (FCA) 
based on nine microsatellite loci in A. schrenckii, A. sinensis and their 
hybrids. The figure shows the distinction between A. schrenckii, A. sinensis 
and A. sinensis^ x A. schrenckii^. 

Additional file 10: Figure S7. Factorial Correspondence Analysis (FCA) 
based on nine microsatellite loci in H. dauricus, A. ruthenus and their 
hybrids. The figure shows the distinction between H. dauricus, A. ruthenus 
and A. ruthenus'^ x H. dauricusS- 

Additional file 11: Figure S8. Factorial Correspondence Analysis (FCA) 
based on nine microsatellite loci in A. baerii, H. Dauricus and their hybrids. 
The figure shows the distinction between A. baerii, H. dauricus and A. 
baerii^ x /-/. dauricusS. 

Additional file 12: Figure S9. Factorial Correspondence Analysis (FCA) 
based on nine microsatellite loci in A. gueldenstaedti, H. dauricus and their 
hybrids. The figure shows the distinction between A. gueldenstaedti, H. 
dauricus and A. gueldenstaedti'^ x H. dauricusS- 
Additional file 13: Figure S3. Assignation of the 466 sturgeons by 
STRUCTURCE analysis based on nine microsatellite loci in 17 sturgeon 
strains. The figure illustrates the existence of different genetic clusters in 
sturgeon breeds revealed the by the analysis of microsatellite genotyping 
data. 



Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

XMZ and WHW contributed equally to this work. XMZ performed data 
analysis and wrote the manuscript. WHW conceived the original ideas. LML 
and XFM helped to perform the analysis and improved the manuscript. JPC 
conceived the study, made substantial contribution to the interpretation of 
the results, and revised the manuscript. All authors read and approved the 
manuscript. 

Acknowledgements 

This research was supported by the Science and Technology Planning 
Project of Guangdong Province, China (201 1 B090300024) and Supported by 
Program for Innovative Research Team of Cross-discipline and cooperation 
from Chinese Academy of Sciences. 

Author details 

^College of Fisheries, Huazhong Agricultural University, 430070 Wuhan, 
China. ^Guangdong Entomological Institute/South China Institute of 
Endangered Animals, 510260 Guangzhou, China. ^Heilongjiang River Fishery 
Research Institute (HRFRI) of Chinese Academy of Fishery Sciences, 150070 
Harbin, China. 

Received: 28 February 2013 Accepted: 21 June 2013 
Published: 28 June 2013 

References 

1. Abbott R, Albach D, Ansell S, Arntzen JW, Baird SJ, Bierne N, Boughman J, 
Brelsford A, Buerkle CA, Buggs R, Butlin RK, Dieckmann U, Eroukhmanoff F, 
Grill A, Cahan SH, Hermansen JS, Hewitt G, Hudson AG, Jiggins C, Jones J, 
Keller B, Marczewski T Mallet J, Martinez-Rodriguez P, Most M, Mullen S, 
Nichols R, Nolte AW, Parisod C, Pfennig K, Rice AM, Ritchie MG, Seifert B, 
Smadja CM, Stelkens R, Szymura JM, Vainola R, Wolf JB, Zinner D: 
Hybridization and speciation. J Evol Biol 2013, 26:229-246. 



2. Jackson DP: In situ hybridization in plants. In Molecular Plant Pathology: 
A Practical Approach. Edited by Bowles DJ, Gurr SJ, McPhereson M. Oxford: 
Oxford University Press; 1991:163-174. 

3. Leary RF, Allendorf FW, Sage GK: Hybridization and introgression between 
introduced and native fish. Am Fish Soc Symp 1995, 15:91-101. 

4. Rieseberg LH: Evolution: replacing genes and traits through 
hybridization. Curr Biol 2009, 19:R1 19-R122. 

5. Wright S: The roles of mutation, inbreeding, crossbreeding and selection 
in evolution. In Proceedings of the Vlth International Congress of Genetics: 
24-30 August 1932; Ithaca 1932, 1:356-366. 

6. Mallet J: Hybridization, ecological races and the nature of species: 
empirical evidence for the ease of speciation. Philos T R Soc B 2008, 
363:2971-2986. 

7. Saetre GP: Hybridization is important in evolution, but is speciation? 

J Evol Biol 2013, 26:256-258. 

8. Soltis PS, Soltis DE: The role of hybridization in plant speciation. 
Annu Rev Plant Biol 2009, 60:561-588. 

9. Abbott RJ, Ritchie MG, Hollingsworth PM: Introduction. Speciation in plants and 
animals: pattern and process. Philos TRSocB 2008, 363:2965-2969. 

10. Swan AA, Kinghorn BP: Evaluation and exploitation of crossbreeding in 
dairy cattle. J Dairy Sci 1 992, 75:624-639. 

1 1 . Mirkena T Duguma G, Haile A, Tibbo M, Okeyo A, Wurzinger M, Solkner J: 
Genetics of adaptation in domestic farm animals: A review. Livest Sci 
2010, 132:1-12. 

12. Herrero-Medrano JM, Megens HJ, Crooijmans RP, Abellaneda JM, Ramis G: 
Farm-by-farm analysis of microsatellite, mtDNA and SNP genotype data 
reveals inbreeding and crossbreeding as threats to the survival of a 
native Spanish pig breed. Anim Genet 2013, 44:259-266. 

13. Stuber CW: Heterosis in plant breeding. Plant Breeding Reviews 1994, 
12:227-251. 

14. Birchler JA, Yao H, Chudalayandi S, Vaiman D, Veitia RA: Heterosis. 
Plant Cell 2010, 22:2105-2112. 

15. Cundiff LV: Experimental results on crossbreeding cattle for beef 
production, j /\n/m Sci 1970, 30:694-705. 

16. Raymakers C: CITES, the Convention on International Trade in 
Endangered Species of Wild Fauna and Flora: its role in the conservation 
of Acipenseriformes. J AppI Ichthyol 2007, 22:53-65. 

17. Wei QW, Zou Y, Li P, Li L: Sturgeon aquaculture in China: progress, 
strategies and prospects assessed on the basis of nation-wide surveys 
(2007-2009). J /\pp/ Ichthyol 201 1, 27:162-168. 

18. Boyer MC, Muhlfeld CC, Allendorf FW: Rainbow trout (Oncorhynchus 
mykiss) invasion and the spread of hybridization with native westslope 
cutthroat trout (Oncorhynchus clarkii lewisi). Can J Fish Aquat Sci 2008, 
65:658-669. 

19. Lajbner Z, Slechtova V, Slechta V, Svatora M, Berrebi P, Kotlik P: Rare and 
asymmetrical hybridization of the endemic Barbus carpathicus with its 
widespread congener Barbus barbus. J Fish Biol 2009, 74:418-436. 

20. Dudu A, Suciu R, Paraschiv M, Georgescu SE, Costache M, Berrebi P: Nuclear 
markers of Danube sturgeons hybridization. IntJMol Sci 201 1, 12:6796-6809. 

21. Schrey AW, Boley R, Heist EJ: Hybridization between pallid sturgeon 
Scaphirhynchus albus and shovelnose sturgeon Scaphirhynchus 
platorynchus. J Fish Biol 201 1, 79:1828-1850. 

22. Ivanova NV, Zemlak TS, Hanner RH, Hebert PDN: Universal primer cocktails 
for fish DNA barcoding. Mol Ecol Notes 2007, 7:544-548. 

23. May B, Krueger CC, Kincaid HL: Genetic variation at microsatellite loci in 
sturgeon: primer sequence homology in Acipenser and Scaphirhynchus. 
Can J Fish Aquat Sci 1 997, 54:1 542-1 547. 

24. McQuown EC, Sloss BL, Sheehan RJ, Rodzen J, Tranah GJ, May B: Microsatellite 
analysis of genetic variation in sturgeon: new primer sequences for 
Scaphirhynchus and Acipenser. TAm Fish Soc 2000, 129:1380-1388. 

25. Esselink GD, Nybom H, Vosman B: Assignment of allelic configuration in 
polyploids using the MAC-PR (microsatellite DNA allele counting-peak 
ratios) method. Theor AppI Genet 2004, 109:402-408. 

26. Pritchard JK, Stephens M, Donnelly P: Inference of population structure 
using multilocus genotype data. Genetics 2000, 155:945-959. 

27. Pastorini J, Zaramody A, Curtis DJ, Nievergelt CM, Mundy Nl: Genetic 
analysis of hybridization and introgression between wild mongoose and 
brown lemurs. BMC Evol Biol 2009, 9:32. 

28. Clark LV, Jasieniuk M: POLYSAT: an R package for polyploid microsatellite 
analysis. Mol Ecol Resour 201 1, 1 1:562-566. 



Zhang et at. Genetics Selection Evolution 2013, 45:21 
http://www.gsejournal.0rg/content/45/l/2i 



Page 10 of 10 



29. Hebert PDN, Cywinska A, Ball SL: Biological identifications through DNA 
barcodes. Proc Biol Sci 2003, 270:31 3-321 . 

30. Hebert PDN, Stoeckle MY, Zemlak TS, Francis CM: Identification of birds 
through DNA barcodes. PLoS Biol 2004, 2:e312. 

31 . Ward RD, Zemlak TS, Innes BH, Last PR, Hebert PDN: DNA barcoding Australia's 
fish species. Philos Trans R Soc Lond B Biol Sci 2005, 360:1 847-1 857. 

32. Costa FO: Jeremy R deWaard, Boutillier J, Ratnasingham S, Dooh RT, 
Hajibabaei M, Hebert PDN: Biological identifications through DNA 
barcodes: the case of the Crustacea. Con J Fish Aquot Sci 2007, 
64:272-295. 

33. Jenneckens I, Meyer JN, Debus L, Pitra C, Ludwig A: Evidence of 
mitochondrial DNA clones of Siberian sturgeon, Acipenser baerii, within 
Russian sturgeon, Acipenser gueldenstaedtii, caught in the River Volga. 
Ecol Lett 2000, 3:503-508. 

34. Hoskin CJ, Higgle M, McDonald KR, Moritz C: Reinforcement drives rapid 
allopatric speciation. Nature 2005, 437:1353-1356. 

35. Vasil'eva ED, Vasil'ev VP, Ponomareva EN, Lapukhin YA: Triple hybrids 
obtained by artificial hybridization of the Russian sturgeon Acipenser 
gueldenstaedtii with the hybrid of the starred sturgeon A. stellatus and 
the great sturgeon A. huso (Acipenseridae): The kind of inheritance of 
some morphological characters and fertility of the parental hybrid form. 
Jlchthyol 2010, 50:605-617. 

36. Verspoor E, Hammart J: Introgressive hybridization in fishes: the 
biochemical evidence. J Fish Biol 2006, 39:309-334. 

37. Richly E, Leister D: NUMTs in sequenced eukaryotic genomes. 
Mo I Biol Evol 2004, 21:1081-1084. 

38. Bensasson D, Zhang DX, HartI DL, Hewitt GM: Mitochondrial pseudogenes: 
evolution's misplaced witnesses. Trends Ecol Evol 2001, 16:314-321. 

39. Antunes A, Ramos MJ: Discovery of a large number of previously 
unrecognized mitochondrial pseudogenes in fish genomes. 
Genomics 2005, 86:708-717. 

40. Tranah G, Campton DE, May B: Genetic evidence for hybridization of 
pallid and shovelnose sturgeon. J Hered 2004, 95:474-480. 

41. Ludwig A, Lippold S, Debus L, Reinartz R: First evidence of hybridization 
between endangered sterlets (Acipenser ruthenus) and exotic Siberian 
sturgeons (Acipenser baerii) in the Danube River. Biol Invasions 2009, 
11:753-760. 

42. Huang D: Species, distributions and conservation of Acipenseriformes. 

Tsinghua Sci Technol 2002, 7:41 6-420. 

43. Campton DE: Natural hybridization and introgression in fishes: methods 
of detection and genetic interpretations. In Population Genetics and Fishery 
Management. Edited by Ryman N, Utter F. Seattle: University Washington 
Press; 1987:161-192. 

44. Doroshov SI: The escape of cultured sturgeon and the interbreeding with 
wild stocks. Proceedings of the Florida Sturgeon Culture Risk Assessment 
Workshop 2000:6-7. 

45. Kelaita MA, Cortes-Ortiz L: Morphological variation of genetically 
confirmed Alouatta Pigrax A. palliata hybrids from a natural hybrid zone 
in Tabasco, Mexico. Am J Phys Anthropol 2013, 150:223-234. 

46. Welsh A, Hill T, Quinlan H, Robinson C, May B: Genetic assessment of lake 
sturgeon population structure in the Laurentian Great Lakes. N Am J Fish 
Mana 2008, 28:572-591. 

47. Notter DR: The importance of genetic diversity in livestock populations 
of the future. J Anim Sci 1 999, 77:61 -69. 



doi:1 0.1 186/1297-9686-45-21 

Cite this article as: Zhang et al.: Genetic variation and relationships of 
seven sturgeon species and ten interspecific hybrids. Genetics Selection 
Evolution 2013 45:21. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



